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ABSTRACT 

An application of the Darboux transformation on a cnoidal wave background in 
the coupled nonlinear Schrodinger equation gives a new solution which describes a 
soliton moving on a cnoidal wave. This is a generalized version of the previously 
known soliton solutions of dark-bright pair. Here a dark soliton resides on a cnoidal 
wave instead of on a constant background. It also exhibits a new types of soliton 
solution in a self-focusing medium, which describes a breakup of a generalized dark- 
bright pair into another generalized dark-bright pair and an "oscillating" soliton. We 
calculate the shift of the crest of the cnoidal wave along a soliton and the moving 
direction of the soliton on a cnoidal wave. 
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1 Introduction 



In this paper, we consider a system of coupled nonlinear Schrodinger (CNLS) equa- 
tions [1, 2, 3] 

d s ip k = -id 2 z ip k -2i<r(\ip 1 \ 2 + \fo\ 2 )ij, k -, k = 1,2. (1) 

These equations, which were named as the Manakov model, are important for a 
number of physical applications like multi- frequency and/or two different polariza- 
tions of light in a fiber [4, 5, 6]. Nowadays, there exist vast amount of exact solutions 
of this system including exact vector solitons [1, 5], bound solitary waves [7], and 
periodic solutions [7, 8, 9, 10]. Especially, quasi-periodic solutions in terms of N- 
phase theta functions for the Manakov model are derived in [11], while a series of 
special solutions are given in [12, 13, 14, 15]. 

Recently, this equation attracts new attention as it can describe the localized 
states in optically induced refractive index gratings [16, 17, 18]. Strong incoherent 
interaction of such a grating with a probe beam facilitates the formation of a noble 
type of a composite optical soliton, where one of the components (described by ip\ 
in Eq. (1)) creates periodic photonic structure, while the other component (^2) 
experiences Bragg reflection from this structure and forms gap solitons. Optically 
induced lattices are very exciting development which open up creating dynamically 
reconfigurable photonic structures. It was shown that the CNLS equation describes 
the propagation of two incoherently interacting beams in a photorefractive crystal 
in the limit of weak saturation regime [19, 18]. 

To analyze the behavior of the soliton in the optically induced lattices, we need 
a solution of the CNLS equation that describes a soliton moving on a cnoidal wave. 
The integrability of the CNLS equation shown by Manakov provides us to obtain 
various type of solutions. Especially, one can obtain the bright soliton in a focusing 
medium by applying the inverse scattering method (ISM) [20]. On the other hand, 
in the case of nonvanishing background fields, the inverse scattering method is tech- 
nically highly involved and only the dark solitons of the simplest one-component 
case have been found in this way. In the two component case, the Hirota method 
has been adopted to obtain dark solitons instead of the ISM [21, 22]. The Hirota 
method, however, does not provide a way to construct a soliton solution moving 
on a cnoidal wave background. Interestingly, the (soliton + cnoidal wave) solution 
can be obtained from the general quasi-periodic solutions of N-phase theta func- 
tions, by taking degenerate limit of the 2-phase solution. In fact, Ref. [23] applied 
this procedure to obtain a solution of the single-component nonlinear Schrodinger 
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equation (NLSE). But solutions from the N-phase theta functions have the so-called 
"effectivization" problem, which is related to extracting the physical solutions by 
taking proper initial conditions [24, 10]. Much more, these solutions have rather 
complicated form, which makes them difficult to apply to real situations. 

In this paper, we employ a simple, but very powerful soliton finding technique 
based on the Darboux transformation (DT). This method is essentially equivalent 
to the ISM, but avoids mathematical technicalities of the ISM. Section 2 introduces 
the main method including the DT and Sym's solution. Section 3 analyzes the 
characteristics of obtained solutions in the case of focusing medium. Section 4 is 
devoted to the case of defocusing medium. Section 5 contains a short discussion and 
Appendices give some proofs of equations in the main text. 



2 The Method 



2.1 Lax Pair 



We first bring the CNLS equation into a matrix form in terms of 3 x 3 matrices E, T 
and E = [T, E], 





E = -aij>* | , T = | -i/2 | , (2) 
such that 

d g E = -d 2 z E + 2E 2 E. (3) 

One can readily check that the components of Eq. (3) are indeed equivalent to the 
Manakov equation in Eq. (1). The signature a is either 1 or —1 depending on 
whether the group velocity dispersion is abnormal (a — 1) or normal (<x = —1), or 
the waveguide is self-focusing (a — 1) or self-defocusing(<r = — 1). One advantage of 
using matrices is that we can write down the associated linear equation (Lax pair); 

(d z + E + AT)^ = 0, {d- z + EE - d z E - XE - A 2 T)^ = 0, (4) 

where A is an arbitrary complex number and ^(z, z, A) is a three component vector. 
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2.2 Darboux Transformation 



The following cnoidal wave solution of the CNLS equation properly describes the 
background periodic wave in optically induced refractive index gratings. For the 
case of focusing medium (a — 1), it is 

r i (z,z)= P Mx,k)e«, ^2 = 0, (5) 

where x — ~p( z ~ vz ), C — [—vz/2 — p 2 (2 — k 2 )z + v 2 z/4\. The defocusing case 
(a = —1) is obtained by substituting p — > ip, k — > ik in Eq. (5) of the focusing case. 
1 Here dn, sn is the standard Jacobi elliptic function. And v is the velocity of the 
cnoidal wave and k e (0, 1) is the modulus of the Jacobi function. As far as elliptic 
functions are involved we employ terminology and notation of Ref. [25] without 
further explanations. To describe the characteristics of the composite optical soliton 
formed on the background grating, we need a soliton solution superposed on the 
cnoidal wave. To obtain a superposed solution of (soliton + cnoidal wave) using 
the Darboux transformation, we need to find a solution of the linear equations (4) 
with ipi, i — 1, 2 given by Eq. (5). We denote the solution of the linear equation as 

(s \ 

a three component vector, \P = si . Then a new solution of (soliton + cnoidal 

\S 2 ) 

wave) is constructed using the Darboux transformation [26, 27, 28], which is 

VP (*, z) = tf(z, z) + <(A - A*) l2j _ aS f l | , 2 , i = 1, 2- (6) 

P0| + C Z^j=l,2 \ s j\ 

Using that Sj satisfy the associated linear equations in Eq. (4), it can be directly 
checked that ipi~ s , i — 1, 2 in Eq. (6) is a new solution of the CNLS equation. 

2.3 Sym's Solution 

Explicitly, the solution of the linear equations in Eq. (4) for a = 1 can be written 
down as 

So = e^ 2 [M^^ 2 (^)^ (^)-iV e -^ 1 (^)^(^)]/^o(^), 



51 = e " </2 [- MelA ^^)^(^) + ^- lA ^(^)^o(^)]/^o(^,, 
s 2 = Ce iXz/2 - ix2 " /2 , (7) 



1 When we substitute p ^ ip only, it becomes a singular solution of the defocusing 
CNLS equation. 



4 



where M, N, C are arbitrary complex numbers and A = —72 + fix with 

p 2 2 , cn 2 (u,kf), 
7 = — ridn (m, fc ) + 



2 v ' ' sn 2 (u,£;') J ' 
d ,m. 1 dn(w, fc')cn(u, fc') sn(u, fc')dn(tx, k') . , 

^ = ^ ln ^ o( 2K ) + 2 s^O + ' () 

and X and are complete elliptic integrals of the first and the second kinds, respec- 
tively. In the above formula, a complex parameter u is related to the DT parameter 
A as following, 

v dn(w, k')cn(u, k') , , 

Sym first introduced this solution in the description of vortex motion in hydrody- 
namics [29]. It was then applied to an NLSE- related problem in Ref. [30]. The proof 
of Sym's solution is given in the Appendix I. Finally, solutions of the case a = — 1 
is obtained by substituting p — > ip, k — > ik in Eq. (7). 



3 Solutions of Self-focusing Medium 
3.1 A Soliton crossing the Cnoidal Wave 

Figure 1 shows (dark soliton + cnoidal wave) |*0i|, and a bright soliton |*0 2 |- It is 
obtained using Eqs. (5), (6), (7), (8) and (9) and taking iV = 0. Other parameters 
for the figure are v = 0, k = 0.9, p = 1.3, u = -0.38 + 0.63i,M = 1,C = 0.3. 
The figure is drawn using MATHEMATICA, which is also used in checking that 
the solution indeed satisfies the equation of motion (1). It shows the characteristic 
dark-bright pair soliton of CNLS equation where the dark soliton resides on a cnoidal 
background. This case becomes the well-known dark-bright soliton pair when we 
take k = [6]. 

The dark-bright pair soliton moves along a line which satisfies | s 2 ! ~ ko^l ~ 
in Eq. (7). Here, |sq M ^| or \s[ M ^\ mean M-coupled terms in Eq. (7). Thus 
Isg^l is obtained from |so| in Eq. (7) by taking M — 1, N — 0. In the subsequent 
section, we use notations \ sq N ^\ or (s^], which means iV-coupled terms in Eq. (7). In 
this case, \sq^\ i s obtained from |so| in Eq. (7) by taking M — 0, N — 1. As we move 
away from the soliton, it becomes \s 2 \ » \sq M ^\ ~ |si M ^| or |s 2 | << l s o M ^| ~ l s i M ^|; 
and ip 2 in Eq. (6) becomes zero. We can see from Eq. (7) that the dominating 
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Figure 1: Dark-bright soliton pair of CNLS equation: j^l shows a dark 
soliton residing on a cnoidal background. \t/j 2 \ shows a bright soliton. The 
parameters are v = 0, k = 0.9, p = 1.3, u = — 0.38 + 0. 63i, M = 1, N = 0, C = 
0.3. 



factor of \sq M ^\ (or Is^l) is exp(— ImA). Other terms give much small oscillating 
behavior. The dominating factor of \s 2 \ is exp(Im(A 2 z — Xz)/2). Thus, the soliton 
pair moves along the line z = az, where 

a = Im(A 2 /2 - 7 + pv(3)/Im((3p + A/2). (10) 

The value a of the soliton for parameters of Figure 1 is calculated to be -1.92 using 
Eq. (10), which is in accordance with Figure 1. 

Another interesting feature of Figure 1 (l^il) is that the crest of the cnoidal 
background shifts constantly across the dark soliton. The shift is calculated as 
following. As we move away from the soliton such that \s 2 \ » |s[) M ^ ~ l s i M ^l; we 
find that ipi~ s (z, z) — ► i\)\{z,z) from Eq. (6). Thus |t/>i~ s | — > pdn% in this region. 
On the other side of the dark soliton, it becomes \s 2 \ << \sq^\ ~ ki^l, and 

^r^,^)-V'i(^^+i(A-A*K4 M) /4 M) +4 M) 74 M) *)- 1 - using that 4 M) /4 M) = 

— e _l ^sn(— iu)dn(x + iu)/cn(— iu) (Use Eqs. (7) and (25) with N = 0), we find that 
(For simplicity, we take u = iuj.) 

,r-s -ic ! c dn-u/ /snw/dn(x - uA cnw 7 \ 1 

ft * e <^pdn X -2p ^ '- + — , 11 

cnui snui \ cnw/ snu/dn^x _ u i) J 

where we use Eqs. (9), (26) and snu/ = sn(-u/, k). Using the addition theorem, Eq. 
(11) can be written 

dn(x — ui)dmii — /c 2 snu/cnu/sn(;\; — w/)cn(x — ui) 



p- 



-2p 



1 — /c 2 sn 2 «/sn 2 (x — ui) 
dn-u/dn(x — uj) 



= P 



sn 2 M/dn (x — Uj) + cn 2 w/ 
dn(x — Ui)dmij — /c 2 snu/cnu/sn(;\; — M/)cn(x — «/) 
1 — A; 2 sn 2 M/sn 2 (x — ui) 
= -pdn(x-2u 7 ). (12) 

Thus, \ipi~ s \ — > pdn(% — 2-u/) in this region, and the shift of crest is 2uj. For complex 
values of u = ur + iui, we find the following identity which is checked numerically 
for various values of complex u and real x, k. 

. / da(-iu,k) dn(iu*,k) \ 

dn fc - t • m < M + c * m V * M x 
\cn(— ra, /c)sn(— ra, k) cn(2-u*, k)sn{tu*, k) J 

/ sn(-m,fc)dn(x + m,fc) | cn(m*,fc) = 

\ cn(— iu, k) sn(iu* , k) dn(x — iu* , k) J 



Using Eq. (13), we find that | Vx ^ I — * £>dn(x-2-u 7 ) = pdn(x-2Imu) still holds. We 
see that the parameter u is suitable (compared to the DT parameter A) in describing 
the characteristics of (soliton + cnoidal) system. The shift of the crest in terms of 
z is z — > z + 2Imu/p when we take v = in x — ~p( z ~ vz)- 111 Figure 1, this shift 
is 0.97. 

3.2 Soliton on top of Cnoidal wave 

Recently there arise new interests on optically-induced lattices and the localization 
of light in those gratings [16, 17, 18]. Especially, Ref. [18] describes a stationary 
configuration of (soliton + cnoidal wave) system where a soliton moves in parallel 
with the crest of a cnoidal wave. In this case, the velocity of the soliton is zero and 
Eq. (10) requires that A 2 /2 — 7 +pvf3 should be a real number. It, in turn, requires 
that u takes forms u = iuj or u = iuj + K' . To see this, use Eqs. (8), (9) and (26). 
Note that [31] 

/ -rw\ 1 / -rw\ i dna . .cna 

sn(a + zif') = , cn(a + = , dn(a + %K') = -1 . (14) 

sna k sna sna 

Here Uj is an real number and K' = K(yl — k 2 ). 

Fig. 2 shows one such case. Here solitons run in parallel with the crest of cnoidal 
wave. For this figure, we take M — 1.1, N — 1. In this case, there arises interference 
between the M and iV-coupled terms in s , si of Eq. (7). This interference results 
in an oscillating behavior of solitons as shown in Fig. 2. 

The oscillating behavior disappears in the case of N = (or M — 0), and we get 
a stationary configuration of (soliton + cnoidal wave) system. This case corresponds 
to the system discussed in [18]. (They use numerical analysis.) Various stationary 
configurations dealt in [18] can be obtained from our result in Eqs. (6) and (7) by 
taking proper values of parameters. It includes cases taking fc>l,0>fc>l,as 
well as taking u = iuj or u = iui + K' with iV = (or M — 0). More detailed 
characteristics of solutions corresponding to this specific choices of parameters are 
discussed elsewhere [32]. 

3.3 Soliton Fusion on a Cnoidal wave background 

More general solution describing (soliton+cnoidal wave) is given by taking M 7^ 
0, iV 7^ in Eq. (7). Figure 3 describes one of this case. We see that a dark-bright 
pair breaks up into another dark-bright pair plus an oscillating soliton. This feature 
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Figure 2: Dark-bright soliton pair moving paralelly along the crest of cnoidal 
wave: for dark soliton, \ip 2 \ for bright soliton. The parameters are v = 
0,A; = 0.9,p= 1.3, u = 0.63i,M = 1.1, N= 1,C = 0.3. 



Figure 3: Dark-bright soliton pair breaks up into another dark-bright pair 
plus an oscillating soliton: shows a dark soliton breaking up into a dark 
plus oscillating solitons. \ip2 | shows a bright soliton changing its direction. 
The parameters are v = 0, k = 0.9, p = 1.3, u = —0.38 + 0.63i, M = 1.1, N = 
1,C = 0.3. 



was first found in [5]. Our result is a generalization of the result in [5], such that a 
dark soliton resides on a cnoidal wave instead on a constant background. 

There are three soliton lines in Fig. 3. The direction of these lines is calculated 
similarly as in the subsection 3.1. The directions of two dark-bright pairs are given 
by lines satisfying either \s 2 \ ~ |sq M ' > | ~ |si M ^| or \s 2 \ ~ l^o^l ~ l s i I- The direction 
of soliton line for \s 2 \ ~ |sq M ^| ~ l s i M ^| is given in Eq. (10) of subsection 3.1. It was 
z = — 1.92z for parameters of Fig. 1 (or Fig. 3). Contrary to the case of Fig. 1 
(where the soliton exists on a full straight line), the soliton in Fig. 3 (which exists 
on a half of the line) exists only in the region satisfying \s 2 \ ~ \ s b M ^\ ~ l s i M ^| >> 
ko^l ~ l s i^l- Using the expression for ip 2 in Eq. (6), it is easy to see that the 
soliton disappears in the region |s 2 | ~ ko^l ~ l s i M ^l << l s o^l ~ l s i^l- 

Another dark-bright pair moves along a line satisfying | «2 1 ' — ' \ s o^\ ~ l s i^|) or 
z = az, where 

Im(AV2 + 7 -^) 

Im(-/?p + A/2) • {0) 

In this case, the dominating factor of \sq N ^\ (or |si^|) in Eq. (7) is exp(ImA). The 
dominating factor of |«2| is still exp(Im(A 2 z — Xz)/2). Thus, the soliton pair moves 
along the line ImA = lm(\ 2 z — Xz)/2), which gives Eq. (15). For parameters of 
Fig. 3, a = 0.44. 

The oscillating soliton in |-^i| of Fig. 3 moves along a line z = az, where 

~ Im(-7+pu/3) 

a = M(3 P ) ■ (16) 

In this case, the line is described by \s M ^\(~ ki^l) ~ l^o^K^ ki^D- The dom- 
inating factor of \sq^\ (or Is^l) is exp(— ImA), while that of \s^\ (or |s^|) is 
exp(ImA). Thus, the oscillating soliton moves along a line ImA = 0, which gives 
Eq. (16). For parameters of Fig. 3, a = —2.89. Note that on a line of the oscillating 
soliton, it becomes \s 2 \ « \s \ ~ |si|. Thus there appears no soliton in \tp 2 \ along 
this line, see Eq. (6). The nature of the oscillating soliton in ip 1 is exactly that of the 
single-component NLSE in [30]. In the region of \s 2 \ » |s | ~ the oscillating 
soliton disappears even in 

The shift of the crest of the cnoidal wave across the solitons can be similarly 
calculated as in the subsection 3.1. Especially in the region (Region I) where \s 2 \ » 
\sq M ^\ ~ |si M ^| and \s 2 \ » \sq N ^\ ~ \s[ N ^\, we have \ipi~ s (z, z)\ — > \ipl(z,z)\ = 
pdn(x). In the region (Region II) where \s 2 \ << \sq M ^\ ~ \s[ M ^\ and \sq^\ ~ 
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\s[ N) \ « |4 M) | ~ \s[ M) \, we have VPM) - + i(A - A*)(4 M) /4 M) + 

5 (m)*^(m)*^_! |-?/>i _s | — > pdn(x — 2iir) = pdn(x — 2lmu), see subsection 3.1. In 
the region (Region III) where \s 2 \ « (4^1 ~ \s[ N ^\ and |4^| ~ (4^1 >> (4^1 ~ 
|4 M) |, we have VPM) -> ViM) + «(A - A*)^/^ + 4^74^*) _1 - In this 
case, Si^/sq^ = — e _l ^cn(— m)/(dn(% — m)sn(— iu)). Similar procedure used in 
Eqs. (11) and (12) gives us \ipi~ s \ — > pdn(x + 2ttj) = pdn(x + 2Inru). 

Region I and Region II meat at a boundary, along which the dark-bright soliton 
pair moves with the velocity a in Eq. (10). Along the boundary between Region 
I and III, the soliton pair moves with the velocity a in Eq. (15). Finally, the 
oscillating soliton moves along the boundary between Region II and III with the 
velocity a in Eq. (16). Thus the relative shift of the crest of the cnoidal wave across 
the boundary described by a is given by x ~^ X ~~ 4Imu. Then the shift of the crest 
in terms of z is z — > z + 4Inru/p when we take v = in x — ~p( z ~ v ^)- For 
parameters of Figure 3, this shift is 1.94. 



4 Solution of defocusing medium 

Figure 4 shows a dark soliton on a cnoidal background in a defocusing medium. 
It is obtained using Eqs. (5), (6), (7) and taking a = —1, v — 0, k — 0.9i,p = 
1.3i, u = -0.3 - 0.98i, M = 1,N = 0,C = 0.1. Note that in the case of defocusing 
medium (er = —1), we take imaginary k,p. Another specific feature of the defocusing 
medium is that it is not permitted to take both M and iV to be simultaneously 
nonzero. This fact is explained as following. The denominator part of Eq. (7) 
is | s 1 2 _ l s i| 2 _ l s 2| 2 ) which is required to be positive definite or negative definite 
on all (z, z) plane. Otherwise, we will get an unphysical singular solution. But 
there always exists a region in the (z,z) plane where |«2| >> |so|, which means 
the denominator must be negative definite. Thus we need |so| 2 — |si| 2 should be 
negative definite. In the region ImA ~ — oo where \sq\ ~ |4 M ^I; ~ 1 4^ 1 5 the 
negative definiteness requires 1 < |4 M V4 M ^I = |sn(— m)dn(% + m)/cn(— m)|. On 
the other hand, in the region ImA ~ oo where |«o| ~ 14^1) l s il ~ 14^1) we need 
1 < (4^/4^1 = l cn ( — ^)/(dn(x — iu)sn(— iu))\. As |dn(%±m)| is an oscillating 
function of x, the negative definiteness from the two regions results in a contradiction 
on the magnitude of |sn(— m)/cn(— iu)\. Thus we need one of M, N should be zero. 
Thus there does not exist the phenomenon of soliton fusion in the case of defocusing 
medium. 
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Figure 4: Dark-bright soliton pair moving in a defocusing medium: 
a dark soliton and 1^2 1 2 , a bright soliton. The parameters are v — 
0.9i,p= 1.3i,u= -0.3-0.98i,M= l,iV = 0,C = 0.1 



When we take N = as in Fig. 4, we need 1 < |sn(— iu)dn(x + iu)/cn{— iu))\. 
In the appendix II, we prove that 



[1 + k 



where 



2\(l/4) 



sn(— iuR + mK', ik) 



cn(—iuR + mK', ik) 



K' 



\dn(x + iu R -mK',ik)\/(l + k 



2\(l/4) 



1, (17) 



(18) 



m is an odd number (m = 2n + 1) 
(Here, we use the fact that k — * 
found that 



K'(l/y/T+k?) 

X = —ip(z — vz) and fc, p, t>, Mr are real numbers. 
ik,p — * ?p in the defocusing medium.) We also 



(1 + k 2 ) 



when 



2\(l/4) 



sn 



-IUr + ttj, i/c) 



cn 



-iur + -u/, i/c) 



>1, |dn( X + iu R - u h ik)\/{\ + k 2 Y 1/4) >1, (19) 



(4n + 1)X' < «/ < (4n + 3)^' (20) 

where n is an arbitrary integer. On the outsides of interval in Eq. (20), Eq. (19) 
reverses its inequality sign. Thus we can see that the negative defmiteness constrains 
the possible value of ui to the intervals given in Eq. (20). Similarly we must take 
the value uj in the outsides of intervals in Eq. (19) when we take M — 0. For the 
parameters of Fig. 4, K' = 0.67. 



The dark-bright pair soliton moves along a line, which satisfies | s 2 1 



I AM) 



(N = case). Thus we can use Eq. (10) to describe the direction of the 
soliton line. In the case of defocusing medium, we must take imaginary p, k values. 
For parameters of Fig. 4, a — —2.48. The soliton line shifts along the z-axis (or 
z-axis) according to the relative magnitude of M and C. This shift is rather large 
in Fig. 4 compared to the other figures, which is due to large M/C = 10 value. 

The shift of the crest for the case of defocusing medium is similarly calculated as 
in subsection 3.1. At the region |s 2 | >> |so^l> ipi~ s {—pi{z — vz)) = ipl(—pi(z — vz)). 



To calculate ipi 

\dn(ix,ik) - 



in the region the \s2\ << 
dn(— iu, ik) 



, we use the following formula 
dn(iu*, —ik) 



^cn(—iu, ik)sn(—iu, ik) cn(iu*, —ik)sn(iu* , —ik) 



x 



'sn(— iu, ik)dn(ix + iu, ik) 

cn(—iu,ik) 
= |dn(i% + 2iReu, ik)\, 



cum 



-ik) 



-i 



sn(iu*, —ik)dn(—ix — iu*, —ik) 



(21) 
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where u is complex and Xi k are real. Then using Eqs. (6) and (21), we find that 
\ipi~ s {— pi{z — vz))\ = \ip\(—pi(z — vz) + 2iR,eu) \ for real p. The shift of the crest in 
terms of z is z — > z — 2Reu/p for the case of v — 0. In Figure 4, this shift is 0.46. 

5 Conclusions 

In this paper, we have introduced (soliton+cnoidal wave) solutions of the CNLS 
equation. They were obtained using the DT and introducing a solution of the asso- 
ciated linear problem on a cnoidal background. We calculate the moving direction 
of a soliton on a cnoidal background and the shift of the crest of a cnoidal wave. We 
also found a solution in the self-focusing case where a dark-bright pair breakup into 
another dark-bright pair and an oscillating soliton. These type of solutions, which 
can be easily applicable to the analysis of physically interesting processes, seem 
rather rare in the literature of physics. These solutions can be used, for example, in 
describing the localized states in optically induced refractive index gratings. 

The stability analysis of these solutions is remained for future study. In fact, 
there already appears some numerical studies on this subject [18]. There was shown 
that the (soliton+ cnoidal wave) system is unstable or weakly stable for the focusing 
case, while it is stable for the defocusing case. 
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6 Appendix I: Proof of Sym's solution 

In this appendix, we show that sq, si, S2 in Eq. (7) indeed satisfy the linear equation 
(4). Consider the following equation, which is obtained from the <9 z -part of Eq. (4); 

d z s + if>° Sl + i\s /2 = 0. (22) 
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By inserting s , Si (For simplicity, we consider a case of M — 1, N — in Eq. (7).) 
into Eq. (22) and taking as in Eq. (5), we get 

_• a V go(Sr) i g ^) 1 , n^IlMS!, dn( M ,fc')cn( M ,fc') 
«PP 2X0 O (^) + 2K0 O (^) P U ' ^ 2 (ff)^o(^) P sn(u,k>) ■ 

(23) 

Inserting (3 in Eq. (8) into Eq. (23) and using following identities [25, 31], 

u 1 Q'(JL.\ E 

s„ U = -)_|!<A> , cn U = ±f dn U = , (25) 

, sn(u,k) / ,a 1 / /x dn(w, fc) 

sn(wt, A; ) = 2 — 7 — — , cn(iu, k ) = — - — — , dahu, k ) = — - — — , (26) 
cn(u, k) cn(u, k) cn(u, k) 

we get 

r\-\-iu sn( — iu k^j 

—p / [dn 2 v — dn 2 (t> — iu)\dv H — , ' [dn(— iu, k) — dn(%, k) dn(x + iu, k)]. 

Jo cn( %u, K ) 

(27) 

Finally, using the following identity (It is a result of the addition theorem of Jacobi's 
elliptic functions.) 

dn 2 (a — b) — dn 2 a = k 2 — -[cna sna dn(a — b) + dna cn(a — b) sn(a — 6)1 

cno 

sn6 d r , , , ... 

[dna dn(a - 6)], (28) 



cn6 da 

we can see that Eq. (27) becomes zero. Other c^-part of the linear equation (4) 
(including the M = 0, N = 1 case of s , s 1 in Eq. (7)) are similarly proved. 
A <9 2 -part of Eq. (4) is (a = 1 case) 

d s s + i\^ c i\ 2 s - id^l s 1 - A^ c si - i\ 2 s /2 = 0. (29) 

By inserting s , s i (For simplicity, we take M — 1, N — in Eq. (7).) into the left 
part of Eq. (29) and taking ty\ as in Eq. (5), we get 

l -[-p\2 - k 2 ) + - i 7 " + i\) + tp 2 dn 2 X - % -\ 2 

„ (~ iu \ ( Z+llL ) 

+ p [ id z dn X + (A + ^)dn X ] ; !f, ; x 2 + l • (30) 
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Using Eq. (22) and the identity (25), Eq. (30) becomes 

+ip 2 dn 2 x — (^u dn% — ip 2 A; 2 snxcnx — Apdnx) — | — , ' ^| dn(% + iu, k). (31) 

With the help of Eqs. (26) and (9), repeated applications of addition theorem on 
Eq. (31) gives zero, which proves Eq. (29). Other <9—part of the linear equation (4) 
(including the M — 0, N — 1 case of s , si in Eq. (7)) are similarly proved. 

7 Appendix II: Proof of Eq. (17) 

We first start with 

+ {m sn(-iu,ik) = + k 

cn(— iu, ik) 

where k = l/y/l + k 2 with k real and u is a complex value. By inserting u = u R +iK' 
in Eq. (18) and using the addition theorem, Eq. (32) becomes 

sn(-Vl + k 2 u R , k)(l + k) - icn(-y/l + k 2 u R , fc)dn(— \/l + k 2 u R , k) 

1 + hn(-y/T+Wu R ,k) ' 

where we use 



sn 
cn 



-iVT+WK', k) = -i(l + A; 2 ) (1/4) , &n(-iVT+k?K', k) = y/l + l, 
(-tVYT¥K', k) = (1 + A; 2 ) (1/4) ^1 + jfc. (34) 



It is now easy to see that the absolute value of Eq. (33) is 1. 
In the case of dn formula, 



dn(-m, ik)/{\ + A; 2 ) (1/4) = cnC-VTTPtt, k) 

dn{-VTT¥u,k){i + k 2 Y 1 /^ 



Applying a similar procedure used in obtaining Eq. (33) to Eq. (35) gives 

cn(— y/l + k 2 u R , k) + isn(— Vl + k 2 u R , fc)dn(— a/1 + k 2 u R , k) 
dn(— y/1 + A; 2 -u fi , fc) — i/ccn(— Vl + fc 2 «ij, k)sn(— y/l + /c 2 w R , fc) ' 
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(36) 



whose absolute value is 1. It is easy to expand above results to the case of u — 
lift ~ iK' . Using that sn(—iu,ik)/cn(—iu,ik),dn(—iu,ik) are periodic functions in 
Ui with the period 2K', we can obtain the result in Eq. (17) in subsection 4. 

To obtain a result shown in Eq. (19), we should go through similar procedures in 
the above paragraphs with much complex expressions in this case. Instead of going 
through a tedious calculation, we use MAPLE to check numerically the validity of 
Eq. (19). 
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